Search for bottleneck effects in Penna ageing and Schulze language model 
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No influence was seen when in two models with memory effects the populations were drastically 
decreased after equilibrium was established, and then allowed to increase again. 
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PACS numbers: 89.65.S, 87.23. Ge, 87.23. Cc 



I. INTRODUCTION 

In reversible models like the Ising paramagnet above 
the Curie temperature, the equilibrium distribution is in- 
dependent of the initial distribution. This is not always 
the case in systems with memory, for example if a small 
fraction of elements always remains at the initial states. 
In biology, at a "bottleneck" (see e.g. [J) most of the 
population dies out e.g. due to an environmental catas- 
trophe, and then the population grows back to about its 
old size. Immediately after this bottleneck, the distribu- 
tion of genes then is an equilibrium distribution, which 
in general differs from both a random distribution and 
the case where all genomes are identical. Then, if the 
system has long-time memory, the final genetic distribu- 
tion long after the bottleneck can be different from the 
one immediately before the bottleneck. 

The present note searches unsuccessfully for such bot- 
tleneck effects in the Penna model for biological ageing 
and the Schulze model for human languages. Though 
both models are reviewed e.g. in ^] we define them in 
Section^land present our simulations in Sect ion UTTI with 
conclusions in Section Hvl 



II. MODELS 

In the Penna model 0, Q , individuals age according 
to a bit-string representing their genome. A bit set to 
zero represents health, a bit set to one represents an in- 
herited disease which starts at the age corresponding to 
the position of this 1 within the 32 bits of the string. 
Three or more active diseases kill the individual. At and 
above the reproductive age of eight (in time units of bit 
positions) each individual has three offspring having the 
same genome except that one randomly selected bit is set 
to one. If it is already one it stays at one. This latter 
assumption makes the model irreversible: Bits set to one 
never can become zero again. Initially, all bits are set to 
zero. 
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The Schulze model also started with bit-strings, but 
later it was generalized to F features each of which 
can be 1, 2, 3, . . . , Q, with F = Q = 5 in our sim- 
ulations below. This string of small integers represents 
the language (or grammar) of the individual. Each site 
of a square lattice carries one individual, which on each 
sweep through the lattice is replaced by a child having 
the same language except for one mutation. And also at 
each such iteration each individual switches to the lan- 
guage of a randomly selected neighbour with probability 
0.9(1 — where x is the fraction of people in the whole 
population speaking the individual's language. The mu- 
tations happen for each of the F features with probability 
p; if a mutation happens, then with probability 1/2 the 
feature gets a random value between 1 and Q, or gets 
the value of a randomly selected nearest neighbour on 
the lattice. These mutations, in contrast to those of the 
Penna ageing model, are reversible. However, the model 
shows a strong first-order phase transition with hystere- 
sis, between dominance of one language and fragmenta- 
tion of the population into numerous languages. Due 
to this hysteresis the results also in this language model 
may depend on the initial conditions or more generally 
on history. 



III. SIMULATION 

In both models we made 10,000 iterations (sweeps 
through the population), then simulated a bottleneck by 
reducing the population to ten randomly selected indi- 
viduals, and then made another 10,000 iterations to find 
the new equilibrium. Then we compare the results after 
the second half with those after the first half (before the 
bottleneck) . 

Fig. n shows for the Penna ageing model that the 
effects of the bottleneck vanish after several dozen iter- 
ations, before even the total population gets its equilib- 
rium value. Also in the age distribution of the population 
and the distribution of 1 bits among the 32 bit positions 
no significant difference was found (not shown.) 

Fig. [3 shows the size distribution of languages, i.e. we 
count how many languages exist which are spoken by a 
given number of people. Again, no effect of the bottle- 
neck is seen. The same holds for the time development of 
the size of the largest language (not shown). (The bot- 
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FIG. 1: Penna ageing model: Total population for short (part a) and long (part b) times, with final agreement before (+) and 
after (x) the bottleneck. 
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FIG. 2: Schulze language model: Histogram of language sizes 
on 1000 X 1000 lattice at mutation probability p = 0.1. The 
lines give data binned by powers of two in the language size 
— number of people speaking one language. 

tleneck consisted in randomly selecting ten survivors and 



to let them randomly populate the whole lattice again.) 



IV. SUMMARY 

These simulations showed no long-time bottleneck ef- 
fects. This does not exclude that with other parameters 
or in other models such effects exist. And for short times 
the bottleneck effects are clearly visible, Fig. ^a). Thus, 
if for example the species Homo Sapiens survived a bot- 
tleneck 10^ years ago, it could be that some genes mutate 
with a probability of less than 1/10^ per year and thus 
we are not yet in genetic equilibrium. 
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